Light S hypernuclei in four-body cluster models 
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Detailed structure calculations in i-Be, glH, ^-Li, =-H and i-Li are performed within the 
framework of the microscopic two-, three- and four-body cluster models using the Gaussian Expan- 
sion Method. We adopted effective HA*' interactions derived from the Nijmegen interaction models, 
which give rise to substantially attractive H-nucleus potentials in accordance with the experimental 
indications. g-H and ^-Li are predicted to have bound states, we propose to observe the bound 
states in future {K~,K'^) experiments using ^Li and ^"B targets in addition to the standard ^^C 
target. The experimental confirmation of these states will provide information on the spin- and 
isospin- averaged HA'' interaction. 



I. INTRODUCTION 

In studies of nuclear interactions, two-body scattering 
data are the primary input for characterizing interaction 
models. However, S = —1 hyperon (Y)-nucleon (N) scat- 
tering data are very limited because of experimental is- 
sues. For 5* = —2 interactions such as AA and EN, there 
are currently no scattering data. Therefore, the existing 
YN and YY interaction models have a substantial de- 
gree of ambiguity. Some YN scattering experiments will 
be performed at the Japan Proton Accelerator Research 
Complex (J-PARC) in the near future. Even at this fa- 
cility, however, the possibility of performing 'E.N or A A 
scattering experiments is very limited or practically im- 
possible. Hence, in order to obtain useful information on 
S = —2 interactions, studies of many-body, hypernuclear 
structure are indispensable. 

Our intention in this work is to investigate the possible 
existence of S hypernuclei and to explore the properties 
of the underlying EN interactions. Identification of E 
hypernuclei in coming experiments at J-PARC will con- 
tribute significantly to understanding nuclear structure 
and interactions in S* = —2 systems, which can lead to 
an entrance into the world of multi-strangeness. In order 
to encourage new experiments seeking E hypernuclei, it 
is essential to make a detailed theoretical investigation 
of the possible existence of bound states, despite some 
uncertainty in contemporary EN interaction models. 

We investigate here the binding energies and struc- 
ture of E hypernuclei produced by {K~ , K~^) reactions 
on light targets on the basis of microscopic cluster mod- 



els. One of the primary issues is how to choose the EN 
interaction. Although there are no definitive data for 
any E hypernucleus at present, a few experimental data 
indicate that S-nucleus interactions are attractive. One 
example is the observed spectrum of the (K~,K~^) re- 
action on a -'^^C target, where the cross sections for E~ 
production in the threshold region can be interpreted by 
assuming a S-nucleus Wood-Saxon (WS) potential with 
a depth of ~ 14 MeV [l|. Other indications of attrac- 
tive S-nuclcus interactions are given by certain emulsion 
data, the events for twin-A hypernuclei, where the initial 
E~ energies were determined by the identification of all 
fragments after the E~p-AA conversion in nuclei. The in- 
ferred E~ binding energies are substantially larger than 
those obtained using only the Coulomb interaction Q]. 
When these E~ states are assumed to be Ip states, the 
WS potentials obtained from the binding energies are 
similar to the one above. These data suggest that the 
average EN interaction should be attractive, which we 
utilize to select the appropriate interaction models. In 
this work we adopt two types of EN interactions, the Ni- 
jmegen Hard-Core model D (ND) ^SJ and the Extended 
Soft-Core model (ESC04) 

The structure of light p-shell nuclei can be reason- 
ably described in terms of cluster models composed 
of two- or three-body subunits. Here, we model the 
possible E~ hypernuclei produced by {K^,K^) reac- 
tions on available light p-shell targets as four-body clus- 
ter structures: The possible targets ^^C, "^^B, ^"B, 
^Be and ^Li naturally lead to such cluster configura- 
tions as aatE~{i^Be), aa2nE~ Q^Li), aanE~ (i^Li), 
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ain5~(-.?.He) and Q;nnS~(_lH), respectively, by conver- 
sion of a proton into a . (In our model calculations, the 
aS~ potential is generated from a G-matrix 'EN inter- 
action via a folding procedure.) Here, among the above 
"E." hypernuclei, -lH(Q;nnS~) is the lightest S~ bound 
system, as shown in the following section. In the case of 
lighter targets, ^Li, ^He, '^He and d, the S~-hypernuclear 
states are composed of anS~, pnnEr pnEr and 

nS~ configurations, respectively. However, these systems 
are not expected to support bound states, considering 
the weakly attractive nature of the EN interactions sug- 
gested so far, except for Coulomb-bound (atomic) states. 
Thus, possible E~ hypernuclear states to be investigated 
lie in the light p-shell region and may be considered to 
have basically a four-body cluster structure. 

This paper is organized as follows: In Sec. II, we de- 
scribe the basic properties of the EN interaction models 
and make clear what is relevant in the present four-body 
calculations. In Sec. Ill, we perform the calculation of 
i^Be(aatS~) with some approximations, in order to fix 
the EN interaction strengths to be consistent with the 
{K~,K^) data. The four-body cluster models, based 
on the Gaussian Expansion Method (GEM), have been 
developed in a series of works for A and double-A hyper- 
nuclei [1, 0, i, i, [13, El, m • In this work, similar cluster 
models are applied to i?Be(aQ;fS~), _lH(Q;nnS~) and 
i°Li(aanS~). In Sec. IV, first we show the calculated 
behavior of the jLH{aE^) and -^Li(aaS^) systems, as 
a function of the kp parameter in the EN G-matrix inter- 
action, to confirm the binding mechanism before adding 
neutron(s). Then, in Sec.V we discuss the calculated re- 
sults for --lH(annS~) and Li(aanS~). 



II. ■EN INTERACTIONS 

As stated above, the experimental information on EN 
interactions is quite uncertain. It should be comple- 
mented by theoretical considerations. Various S'f/s-based 
interaction models have been proposed so far. In the 
construction of these models, the scarce YN scattering 
data are supplemented by the rich NN scattering data 
through use of SU^ relations among the meson-baryon 
coupling constants. Though these models are more or 
less similar in 5 = — 1 systems, their S = —2 EN pre- 
dictions differ dramatically from one another; most are 
repulsive on average. In order to generate an attractive 
EN interaction on the basis of OBE modeling, it seems 
to be necessary that specific features are imposed. In the 
past, the ND model has been popular for S = —2 inter- 
actions, because this model is compatible with the strong 
AA attraction indicated by the older data on double A 
hypernuclei, and also it yields attractive S-nucleus inter- 
actions. These aspects of ND are the result of its specific 
feature that the unitary-singlet scalar meson is included 
without any scalar-octet mesons. In this case, the strong 
EN attraction originates from this scalar-singlet meson 
which gives the same contributions in all YN and YY 



channels. In the case of other Nijmegen OBE models, 
the attractive contributions of the scalar-singlet mesons 
are substantially cancelled by those of the scalar-octet 
mesons, and their EN sectors are repulsive on average. 
A different OBE modeling for attractive EN interactions 
has been adopted in the Ehime model where the 
insufficient EN attraction given by scalar-nonet mesons 
is supplemented by adding another scalar-singlet meson 
a and the coupling constant g-ssa is adjusted so as to 
give reasonable EN attraction, independent of the SU3- 
relations among coupling constants. The two models, 
ND and Ehime, are essentially similar, in that substantial 
parts of the EN attraction result from the scalar-singlet 
mesons. 

More recently, new interaction models ESC04 (a,b,c,d) 
have been introduced, models in which two-meson and 
meson-pair exchanges are taken into account, and in prin- 
ciple no ad hoc effective boson-exchange potentials are 
included 0,11. The features of the ESC04 models differ 
significantly from those of the OBE models, especially 
in the S = —2 channels. Among the ESC04 models, 
ESC04d is distinguished, because the resulting S-nucleus 
interaction gives attraction suggested by the above ex- 
perimental situation. This is mainly due to the following 
mechanism: A remarkably strong attraction appears in 
the T = triplet-even (^^S":) state, because the strongly 
repulsive contribution of vector mesons is cancelled by 
the attractive contributions from axial- vector mesons. In 
fact, the attraction in this state is so strong that pecu- 
liar E bound states are produced in few-body systems 
[H, though such considerations lie outside the scope of 
the work presented in this paper. In later calculations, 
the important points are the spin- and isospin-averaged 
even-state interactions, which are strongly attractive ow- 
ing to the significant ^^S'l-state attraction. Another im- 
portant feature of ESC04d in the S = —2 channel is 
that the meson-pair exchange terms give rise to strong 
AA-SA^-EE and EN-AT.-T.il couphng interactions. This 
feature of ESC04d makes the conversion widths of E- 
hypernuclear states far larger than those for ND. 

TABLE I: Partial wave contributions to Us.[po)- In the case of 
ESC04d, the medium-induced repulsion is included by taking 
av = 0.18. In the case of ND, the hard-core radii are taken 
as Tc ~ 0.52 and 0.45 fm in the ^^^o and the other states, 
respectively. 



model 


T 


'So 




^Pi 




^Pi 


^P2 


Us. 


Th 


ESC04d 





6.3 


-18.4 


1.2 


1.5 


-1.3 


-1.9 








1 


7.2 


-1.7 


-0.8 


-0.5 


-1.2 


-2.5 


-12.1 


12.7 


ND 





-3.0 


-0.5 


-2.1 


-0.2 


-0.7 


-1.9 








1 


-4.1 


-4.2 


-3.0 


-0.0 


-3.1 


-6.5 


-29.5 


0.8 



Our cluster models are composed of cluster units (a 
and i), n and S~, where the a{t)E~ interactions are 
obtained by folding the EN G-matrix interactions into 
the density of a{t). According to the method described 
in Refs. [J, Q , the EN G-matrix interactions are derived 
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from ESC04d and ND in nuclear matter, where the imag- 
inary parts arise from the energy-conserving transitions 
from SiV to AA channels in the nuclear medium. The 
resulting complex G-matrix interactions are represented 
as fci?-dependent local potentials 

3 

Ci^jir, kp) = {a, + bikp + c^kl) exp (-rV/3f ) , 

i=l 

(2.1) 

where kp is, the Fermi momentum of nuclear matter. The 
suffixes and (— ) specify even and odd, respectively. 
In our applications to finite S systems, it is plausible to 
obtain the kp values from the average density in the re- 
spective systems. In a similar G-matrix approach to A 
hypernuclei, for instance, the AA^ G-matrix interactions 
can be adopted to reproduce the observed A binding en- 
ergy (i?A) by choosing appropriate kp values. Such a 
procedure cannot be applied strictly in the case of S hy- 
pernuclei, because there exist no definitive experimen- 
tal data. In this work, we are obliged to choose the kp 
values rather arbitrarily but within a reasonable range 
(0.8 ~ 1.2 fm~^ in light p-shell systems). Here, the ex- 
perimental indication for the existence of i?Be is used 
to adopt the 'E.N G-matrix interactions, although it is 
not so definite because an experimental S binding en- 
ergy (Bh) could not be extracted. As shown later, the 
adjustable parts included in ND and ESC04d are deter- 
mined so that the S s-state energy in our model of Be 
has the value —2.2 MeV for an adequate value of fcp, 
being obtained from the S-nucleus WS potential with 
the depth -14 MeV [l| where the Coulomb interaction 
is switched off. In the case of ND, this constraint can 
be realized by choosing the hard-core radius Tc'- We take 
Tc — 0.52 and 0.45 fm in the ^^Sq state and the other 
states, respectively. The former choice is made so that 
the derived ^^^o AA G-matrix interaction reproduces the 
A A bond energy observed in the double- A hypernucleus. 
On the other hand, the constraint in the case of ESC04d 
is enforced by changing the parameter av controlling the 
medium-induced repulsion We take ay — 0.18. Here- 
after, ESC04d with ay = 0.18 is denoted as ESC for 
simplicity. 

In the Table HI we show the partial- wave contributions 
of the resulting S potential depth Us. in nuclear matter 
at normal density po (kp = 1.35 fm~^). The C/h values 
are found to be very different for ESC and ND, because 
the odd-state contributions in the former are far more at- 
tractive than those in the latter. It is noted, however, the 
odd-state interactions play minor roles in light systems 
considered in this work. More important is that the spin- 
and isospin-dependence differs significantly between ESC 
and ND. 

The interaction parameters (a^, bi and Ci) in our G- 
matrix interactions (|2.1[) are tabulated in Tables [Til and 
mil for ESC and ND, respectively. Hereafter, G-matrix 
interactions derived from ESC and ND are denoted as 
Gpsc and Gnd, respectively. 

The features of our G-matrix interactions can be 
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FIG. 1: The volume integrals of = (G^J' + 3G[,|' + 

are drawn as a function of kp by solid 
curves in (a) for ESC and in (b) for ND. Here, T = 
(dashed) and T — 1 (dotted) parts show the volume inte- 
grals of (G{,J^ +3G^|V4 and (G^^^' + 3G5|>)/4, respectively. 
In (c), the volume integrals of "^^Si and ^^So components for 
ESC are drawn by solid and dashed curves, respectively. The 
corresponding ones for ND are in (d). 



demonstrated clearly by the volume integrals of the G- 
matrix interaction: Jv{kp) = G{r, kp)r^dr. Here, 
we define the spin- and isospin-averaged interactions as 
G(±) = (g[,o^ + 3G[,i^ + 3G[q '' + 9G^i Vl6- The vol- 
ume integrals of G^'^^r, kp) are drawn as a function of 
kp in Fig. [l] where (a) and (b) are for ESC and ND, 
respectively. It should be noted here that the even-state 
interaction of ESC is more attractive than that of ND. 
In the cases of our cluster systems, S-states are deter- 
mined dominantly by a S folding interactions derived 
from G*^*^(r, kp). The G*^~^ for ND is far more attractive 
than that for ESC, though their contributions in s-shell 
systems are very small. Similarly, (c) and (d) in Fig. [1] 
show the volume integrals of the triplet- and singlet-even 
state interactions in the T = 1 state for ESC and ND, 
respectively. Here, the ^^Si and ^^Sq interactions in ESC 
are found to be attractive and repulsive, respectively. On 
the other hand, both of '^^S'l and '^^Sq interactions are at- 
tractive in ND, and the latter is more attractive than the 
former. Namely, the T = 1 spin-spin interaction in ND 
(ESC) is repulsive (attractive). This difference of the 
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T = 1 spin-spin interactions for ESC and ND is reflected 
in the level structures of _IH and i^^Li, as shown later. 

Another important difference between ESC and ND is 
that the AA-SA^-SS coupling interaction in the former is 
far stronger than that in the latter. This is reflected by 
the fact that the calculated value of the conversion width 
Fh for ESC is far larger than that for ND, as exhibited 
in Table U 

TABLE II: The parameters in the G-matrix interaction 
G^'(r-, kp) given by (I2.1|l for ESC. Entries are given in units 
of a [MeV], b [MeV-fm] and c [MeV-fm^] 





ft (fm) 


0.50 


0.90 


2.00 


r'(+) 

'-'00 


a 
b 

c 


0.0 
0.0 
0.0 


-690.8-309.0j 
1263.+252.4i 
-451.7-lll.Oj 


-2.759 
0.0 
0.0 


'-'01 


a 
b 

c 


-6959. 
11280. 
-4371. 


756.5 
-1567. 
627.2 


-1.317 
0.0 
0.0 


^(-) 

^00 


a 
b 

c 


-1634. 
3426. 
-965.8 


257.8 
-137.4 
60.78 


-1.528 
0.0 
0.0 


^(-) 


a 
b 
c 


-5692. 
7697. 
-2667. 


175.0-15.50i 
-583.9-f24.31i 
303.8-13.91i 


-1.411 
0.0 
0.0 


r<(+) 

'-'10 


a 
b 

c 


-216.4 
676.0 
-198.1 


48.96 
-83.76 
43.36 


-1.838 
0.0 
0.0 




a 
b 

c 


527.9 
85.16 
13.25 


-121.8 
-10.83 
9.351 


-1.787 
0.0 
0.0 


'-'10 


a 
b 

c 


-2671. 
3343. 
-1034. 


36.08 
-116.6 
53.97 


-1.043 
0.0 
0.0 




a 
b 

c 


1435. 
451.1 
-131.2 


-166.3 
-24.19 
18.44 


-1.168 
0.0 
0.0 



Our cluster models for A = 7 and 10 systems are com- 
posed of annEr and aanEr , respectively, where the 
'E.N G-matrix interactions are used to obtain qEt fold- 
ing potentials based on the {Osi/2)^ configuration with 
hj^ — 1.358 fm. It is problematic, on the other hand, 
to use the G-matrix interactions for the En parts. The 
reason is as follows: Correlations of En pairs are treated 
exactly in our model space spanned by Gaussian func- 
tions, which means some double counting for En short- 
range correlations that has been already included in the 
G-matrix interactions. Though a reasonable way out of 
this problem is to use directly the bare potentials (ESC 
and ND), there appear some difficulties in such treat- 
ments: In the case of ESC, the SA^-AE and SA^-AE- 
SE coupling potentials in the T = \ channels make our 



TABLE III: The parameters in the G-matrix interaction 
Cij^J {r, kp) given by (|2.ip for ND. Entries are given in units 
of a [MeV], 6 [MeV-fm] and c [MeV-fm^] 





P^ (fm) 


0.50 


0.90 


2.00 


^(-1-) 

'-'00 


a 
b 

c 


-8769. 
15530. 
-6383. 


456.1- 102. li 
-1082.-f91.45i 

473.2- 18.03i 


-2.505 
0.0 
0.0 


'-'01 


a 
b 

c 


452.4 
-25.82 
67.40 


-105.8 
10.75 
10.22 


-.6861 
0.0 
0.0 


r'(-) 
'-'00 


a 
b 

c 


-7382. 
8672. 
-3145. 


-168.9 
-140.1 
69.22 


-3.141 
0.0 
0.0 


M-) 
'-'01 


a 
b 
c 


-569.3 
2072. 
-696.8 


-231.l-7.788i 
-32.47+6. 124i 
22.21-.5631i 


.0300 
0.0 
0.0 


'-'10 


a 
b 

c 


356.9 
110.3 
-1.818 


-138.5 
13.97 
7.792 


-.3949 
0.0 
0.0 


'^11 


a 
b 

c 


436.0 
6.513 
46.10 


-108.4 
10.11 
10.88 


-1.334 
0.0 
0.0 


'-'10 


a 
b 

c 


75.12 
939.6 
-269.5 


-254.3 
-.0260 
7.792 


.1086 
0.0 
0.0 


'-'11 


a 
b 

c 


-281.0 
1227. 
-422.7 


-218.5 
-5.773 
11.82 


-1.003 
0.0 
0.0 



treatment extremely complicated. In the case of ND, 
although these coupling potentials are not taken into ac- 
count, the hard-core singularities cannot be treated in 
our Gaussian model space. Thus, we adopt here sim- 
ple three-range Gaussian substitutes simulating the bare 
potentials. They are fitted so that the G-matrices de- 
rived from them simulate the original T — I G-matrices 
at kp = 1.0 fm^i. Here, the EN-AT, and SiV-AE-EE 
couplings in the ESC case are effectively renormalized 
into the EN single-channel potentials. The determined 
interaction parameters are given in Table HVl for ESC and 
ND. 



III. gi!Be(aafH-) SYSTEM 

Let us start from the analysis for the -L^Be (^^B-f S^) 
hypernucleus produced by the ^^C(Ar^,Ar+) reaction, 
adopting the aatE'^ four-body model. In this case, the 
{T,J'^) = (1, 1~) states are produced, because the 
component is transformed by ATj = 1 on the T = tar- 
get. This system is important in a double sense. One is 
that in BNL-E885 a fairly deep "B-S" potential was m- 
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TABLE IV: Parameters of the three-range Gaussian interac- 
tions simulating (a) ESC and (b) ND in the T — 1 Hn states. 

(a) ESC 



ft (fm) 


0.40 


0.80 


1.50 


31^ 


519.5 


66.27 


-7.230 


33^ 


217.4 


-170.0 


-7.058 


3lQ 


0.0 


-39.56 


-5.178 


33 o 


0.0 


-55.40 


-6.936 


(b) ND 


ft (fm) 


0.50 


0.90 


2.00 


31^ 


1076. 


-159.6 


-5.432 


33^ 


1331. 


-134.0 


-7.610 


3lQ 


0.0 


-32.30 


-5.432 


33 Q 


0.0 


18.16 


-7.610 



dicated as mentioned in the previous section. The other 
is that this reaction is planned as the Day-1 experiment 
at J-PARC. 

The above-mentioned G-matrix interactions Gesc ^-nd 
Gnd are adjusted so as to be consistent with the Woods- 
Saxon potential depth of BNL-E885 within the frame- 
work of the aatEr four-body model. 



Model and Interaction 





C=1 



C=2 



FIG. 2: Jacobian coordinates for the aafS. (iiBe) four-body 
system. The two a clusters are to be symmetrized. 



In the case of an aatEr four-body model, we take two 
sets of Jacobian coordinates as shown in Fig. [21 since we 
get sufficiently converged energies using only those two 
sets of Jacobian coordinates. The total Hamiltonian and 
the Schrodinger equation are given by 



(i/-£;)*jAf(H-Be) = , 



= r + ^ Kfc + Fpauli 
a, 6 



(3.1) 



(3.2) 



where T is the kinetic-energy operator and Vab is the 
interaction between constituent particles a and h. The 
OCM projection operator VpauU will be given below. The 
total wavefunction is described as a sum of amplitudes 
of the rearrangement channels (c = 1 and 2) of Fig. [2] in 
the LS coupling scheme: 



c=ln,N,vl,L,XS,I,K 

X 5J$(ai)$(a2)[$i(t)xi(2-)] 



nlNLvXSIK 



[[0if(r,)^]^i(R,)],e^p, 



cl\ K 



JM 



X [m(i)'7i(S-)] 



(3.3) 



Here the operator Sa stands for the symmetrization op- 
erator for exchange of two a clusters, xi (2~)is the spin 
function of the 2~ particle and ry^(S~) is the isospin 
function of the particle. Following the Gaussian Ex- 
pansion Method (GEM) [H, [H, [i3| , we take the func- 
tional form of (?!)„im(r), V7VLA/(R) and ^l^l^iPc) as 

V'wlm(R) = i?^e-(^/«")ViM(R) , 

U^.{P) = /e-^^/^-^VA^O), (3.4) 

where the Gaussian range parameters are chosen accord- 
ing to geometrical progressions: 



rn 




{n = 


Rn 




{N = 









iV„,ax) 
^max) • 



(3.5) 



The eigenenergy E in Eq. (|3.ip and the coefficients C 
in Eq.( (|3.3p ) are to be determined by the Rayleigh-Ritz 
variational method. 

As for the aa and at interactions, we employ the po- 
tentials which have been used often in the OCM-based 
cluster-model study of light nuclei: Our potentials Yaa 
and Vat [ill reproduce reasonably well the low- lying 
bound states and low-energy scattering phase shifts of 
the aa and at systems, respectively. The Coulomb po- 
tentials are constructed by folding the p— p Coulomb force 
into the proton densities of all the participating clusters. 

The Pauli principle between nucleons belonging to a 
and a;(= a, <) clusters is taken into account by the orthog- 
onality condition model (OCM) [13]. The OCM projec- 
tion operator VpauU appearing in Eq. (j3.2p is represented 
by 

^Pauii - hm 7 ^ |</./(r„,))(0;(r:,J| , (3.6) 

/ 

which rules out the amplitude of the Pauli-forbidden 
a — X relative states f^fiyax) from the four-body to- 
tal wavefunction [21}. The forbidden states are / = 
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OS', IS, OP, OD for X = t and / = OS, IS, OD for x = a. 
The Gaussian range parameter b of the single-particle Os 
orbit in the a particle (Os)'* is taken to be & = 1.358 
fm so as to reproduce the size of the a particle. For 
simplicity the same size is assumed for the t cluster in 
treating the Pauli principle. In the actual calculations, 
the strength 7 for Vpauii is taken to be 10^ MeV, which 
is large enough to push the unphysical forbidden state to 
the very high energy region, while keeping the physical 
states unchanged. 

Using the Vaa and Vat potentials, we perform the 
three-body calculation for the ^^B{aat) system. The 
calculated values of the ground {5/2^) and the first ex- 
cited (l/2j^) states in ^^B are overbound in comparison 
with the experimental values. In order to put the sub- 
sequent four-body calculations for i?Be(aQ;tS^) on a 
sound basis, we introduce a phenomenological aat three- 
body force of the following form: 

Vaat = voexp[-(r„,t//3)2 - (r^.t//?)^] . (3.7) 

Here we adopt Vq — -l-95MeV and /3 — 2.26 fm in or- 
der to reproduce the *^B(3/2j~) ground state energy. For 
the excitation energy of the ^*B(l/2j~) state, we stick 
to the exact experimental value instead of the calculated 
value, when we perform the hypernuclear four-body cal- 
culations. 



B. Results for i-Be and the appropriate fcp 
parameter 

As mentioned before, our 'E.N interactions are adjusted 
so as to give the S~ s-state energy —2.2 MeV in the 

Be system. This value is consistent with the observed 
spectrum of the ^^C{K^ , K^) reaction which suggests 
the WS potential depth of 14 MeV If we assume 
spin-nonflip dominance for the ^^C(K~ , K~^) reaction, 

the [p~y2 ^i/2\j='L~ state is naturally excited. Therefore, 
within the framework of the aatS^ four-body model, the 
fcp parameters in the aS~ and iS^ potentials, without 
Coulomb interaction, are tuned so that the 1^ state en- 
ergies agree with —2.2 MeV. We listed in Table El the 
calculated 2~-binding energies (-83-) of the 1]~ and 2^ 
states. In the case of Gesc, the 2^ state is obtained at 
a lower energy than the 1]~ state. On the other hand, 
the use of Gnd leads to the opposite order. In our 
model, this is because "^^Si interaction for ESC(ND) is 
more(less) attractive than the ^^Sq interaction as shown 
in Fig. 1(c) and (d). The contribution of the E~a and E^t 
Coulomb forces amounts to about 1.5 MeV. The conver- 
sion widths obtained from the imaginary part of Gesc 
is far larger than that for Gnd- This is because the ^Sq 
AA-SA^-EE coupling interaction in ESC is far stronger 
than that in ND. 

We found the appropriate fcp parameter values of the 
effective EN interactions to be kp = 1.055 fm~^ (ESC) 
and k-p = 1.025 fm~-'^(ND), which are consistent with 



the experimental indication in Be. These interactions 
provide our basis to investigate the A = 7 and 10 E~ 
hypernuclei. 



IV. RESULTS FOR TYPICAL SYSTEMS 
COMPOSED OF aH (g!iH) AND aaH (^^Li) 

Let us study the aS^ and aaE~ systems in order to 
demonstrate the basic features of the aE~ interactions. 
In the cases of -IlI(annS~) and i!?Li(aan2~), the dom- 
inant parts of the E~ binding energies are given by the 
aE~ interactions because of the weak binding of the addi- 
tional neutrons. The aE~ interaction is derived by fold- 
ing the EN G-matrix interaction into the wave function 
of the a. The spin- and isospin-dependcnt parts, being 
remarkably different between ESC and ND, vanish in a 
folding procedure involving a spin- and isospin-saturated 
system such as the a. Thus, the aE~ interaction is deter- 
mined only by the spin- and isospin-averaged EN inter- 
action G^^-* (r; fcp'), where the contribution of odd-state 
part G^~^ is quite small in the two-body aS~ system. It 
should be stressed that a-cluster systems such as aE~ 
and aaE~ give the most basic information on the spin- 
and isospin-averaged parts of EN interactions. These 
parts correspond to so-called spin-independent parts in 
interactions represented by the (era), (rr) and {aa){TT) 
operators. 

Here, it is of vital importance how one chooses the kp 
parameters in our G-matrix interactions. The parameter 
kp specifies the nuclear matter density in which the G- 
matrix interactions are constructed. It is most plausible 
that a corresponding value in a finite system is obtained 
from an average density. Our basic interactions (ESC and 
ND) are adjusted so that the derived G-matrix interac- 
tions give rise to reasonable E~ binding in an A ~ 12 
system for kp — 1.0 ~ 1.1 fm~^ adequately chosen. Con- 
sidering that the suitable values kp = 1.055 fm~*(ESC) 
and 1.025 fm^^(ND) for i^Be, it is a modest change to 
take kp — 0.9 fm^* in the A = 4 ^ 6 systems. In fact, 
we have had successful prior experience. In Ref.Q, we 
studied the structure of ^He, ^Be and using the AA^ 
G-matrix interactions, where consistent results were ob- 
tained by choosing the kp parameters to be around 0.9 
fm ~* for ^He and ^Be and to be around 1.1 fm~^ for 
]^C. Here, we take three values of kp parameters for our 
EN G-matrix interactions in order to study A = 7 and 
A= 10 systems: fcF=0.9, 1.055 and 1.3 fm-* for ESC, 
and A:f=0.9, 1.025 and 1.3 fm"! for ND. So, the kp val- 
ues for A — 10 system are considered to be kp ^ 1.0 
fm~^, while those for A — 6 near kp = 0.9 fm^*. The 
unreasonably large value of kp = 1.3 fm~*, as a trial, 
is used only to demonstrate the kp dependences of the 
results. 

In Table rvTl we show the calculated energies and r.m.s 
radii for the aE~ system for three kp values of the EN 
G-matrix interactions. Of course. Coulomb bound {E~- 
atomic) states are obtained, even if the strong inter- 
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TABLE V: The calculated Bg- (MeV) of the Ij and 2j states using ESC and ND potentials for g-Be. In order to reproduce 
the 'observed' -B=- , we tuned fcp — 1.055(fm~^) and 1.025(fm~^) for the ESC and ND potentials, respectively. The energies 
using the ESC and ND potentials without and with Coulomb potentials between a and S~ and between triton and H~, are 
listed respectively. 



ESC ND 

(without Coulomb) (with Coulomb) (without Coulomb) (with Coulomb) 
1" BH-(MeV) 2.24 4.98 2.23 4.82 

r (MeV) 3^95 4^64 L38 1.66 

2- (MeV) 3.18 6.08 1.56 4.06 

r (MeV) 4.24 4.80 0.93 1.18 



TABLE VL The calculated energies of the 1/2+ state, E, 
and r.m.s radii, r^_^-, in the aH~(glH) system for several 
values of fcp. The values in parentheses are energies when the 
imaginary part of the aE~ interactions are switched off. The 
energies are measured from the a + H~ threshold. 

(a) aH-(ESC) 



with fcF(fm-i) 0.9 1.055 1.3 

Coulomb E (MeV) -1.36 -0.26 -0.14 
(-1.71) (-0.57) (-0.19) 

r(MeV) 2.64 0.86 0.15 

r^-5-(fm) 3.89 6.83 13.6 

without E (MeV) unbound unbound unbound 

Coulomb (-0.64) 

r(MeV) - - 



(b) aH-(ND) 



with itF(fm-i) 0.9 1.025 1.3 

Coulomb E (MeV) -0.57 -0.32 -0.15 
(-0.57) (-0.32) (-0.16) 
r(MeV) 0.16 0.06 0.004 

r^_s-(fm) 6.87 9.82 15.65 

without E (MeV) unbound unbound unbound 
Coulomb r(MeV) - - 



actions arc switched off. If a bound state is obtained 
without the Coulomb interaction, this state is called a 
nuclear-bound state. When a nuclear-unbound state be- 
comes bound with help of the the attractive Coulomb in- 
teraction, such a state is called a Coulomb-assisted bound 
state. Table IVll summarizes our results that in each case 
the lowest state is found to be a Coulomb-assisted bound 
state, namely there appears no nuclear-bound state. It 
is noted, in this Table, that a nuclear-bound state is ob- 
tained in the case of GEScikp = 0.9) if its imaginary 
part is switched off: Though the real part of the aS" 
interaction for GESc{kp = 0.9) is attractive enough to 



give a nuclear-bound state, the strong imaginary part 
makes the resulting state nuclear-unbound. In any case, 
the spin- and isospin-averaged even-state part G(+) for 
Gesc is far more attractive than the corresponding part 
of Gnd- 

In Table IVIIl we list the calculated results for the 
aaS~ system. It should be noted, here, that nuclear- 
bound states are obtained in both cases of Gesc f^nd 
Gnd unless an unreasonably large value of A;f is chosen. 
The calculated energies for Gesc are naturally larger 
than those for Gnd- In the aaS~ system, however, G*^~-* 
contributes significantly. It is remarked that the odd- 
state interaction in Gn d is far more attractive than that 
in GesCi which works to reduce the difference between 
both potentials in the aaS" system. 

One notices in the Tables, that the decay widths for 
Gesc are much larger than those for Gnd, when they are 
calculated for the same value oi kp- This is because the 
imaginary part of Gesc is stronger than that of Gnd- 
The difference of the imaginary parts originates mainly 
from the different strengths of ^"^So AA-SA^— EE coupling 
interactions in ESC and ND. 

As mentioned before, noting that the choice of the 
kp value ~ 0.9 fm~^(Q;S~) and ^ 1.0 fm~^ (aaS~), 
are reasonable, respectively, we can expect the exis- 
tence of nuclear-bound states, especially, in the latter 
case. Thus, we can say that observations of oEr and 
aoEr systems certainly provide information about spin- 
independent parts of the 'E.N interactions. In reality, 
however, there are no corresponding nuclear targets to 
produce the above systems by the {K~,K^) reaction. 
As their actual substitutes, in the following, we investi- 
gate the structures of -7_'i{{ann'E.~) and l^'L\{aan'E.^) 
having additional neutron(s), and propose to perform 
the ^Li(i4r~, K^) and ^°B(isr^, K'^) reaction experiments 
with available targets. 



V. A = 7 AND .4 = 10 H HYPERNUCLEI 

Here, we study -IH and i°Li on the basis of annS~ 
and aanEr four-body cluster models, respectively. In 
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TABLE VII: The calculated energies of the 1/2+ state, E, in 
the aoST {JL\a) system for several values of fcp. The values 
in parentheses are energies when the imaginary part of the 
QfH" interaction are switched off. The energies are measured 
from the q:qH~ three-body breakup threshold. 

(a) aaS.- (ESC) 



with fcF(fm ^) 




0.9 


1.055 


1.3 


Coulomb E (MeV) 




-4.81 


-2.23 


-0.83 






-5.17) ( 


-2.57) 


(-1.04) 


r(MeV) 




5.01 


2.89 


1.18 


without E (MeV) 




-2.54 


-0.41 


unbound 


Coulomb 


( 


-2.94) ( 


-0.77) 




r(MeV) 




4.48 


2.18 




(b 




(ND) 




with fcF(fm ^) 




0.9 


1.025 


1.3 


Coulomb E (MeV) 




-2.87 


-1.82 


-0.79 






-2.89) ( 


-1.83) 


(-0.79) 


r(MeV) 




0.58 


0.3 


0.06 


without E (MeV) 




-1.02 


-0.25 


unbound 


Coulomb 


( 


-1.03) ( 


-0.26) 




r(MeV) 




0.45 


0.20 





cluster-model studies, it is essential that interactions 
among cluster subunits be given consistently with respect 
to the corresponding threshold energies. Namely, low- 
energy bound-state energies and scattering phase shifts 
of an, aa, ,ann, and aan subsystems should be repro- 
duced reasonably by the corresponding interactions. We 
emphasize that these severe constraints are correctly sat- 
isfied in the present models, as mentioned below. 





C=2 



h 

C=3 




C=5 - 




9 © 9 




C=7 



C=8 



C=9 



FIG. 3: Jacobian coordinates for all the rearrangement chan- 
nels (c = 1 ~ 9) of the aanSr {^\A) four-body system. 
Two a clusters are to be symmetrized. In the case of the 
annS~(g.lH) four-body system, the two a clusters are re- 
placed by two neutrons, and the neutron is replaced by an a 
cluster. 



in the LS coupling scheme: 
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c=l n,N,u l,L,X S,T.J,K 

X ^jv[$(a)[xi(S-)[xi(ni)xi(n2)] 

X [['^if(rc)^^l(Rc)],eA'(Pc)]^ 
X [77i(S-)[77i(ni)?7i(n2)] Jj,j,^ , 



\SY.IK 



Si s 



JM 



(5.3) 



A. Model and Interactions 

For -IH and i'^Li, all nine sets of the Jacobian co- 
ordinate of the four-body systems are shown in Fig. [31 
respectively. The total Hamiltonian and the Schrodinger 
equation are given by 



[H -E)^.jm{J-^, H°Li) =0 



(5.1) 
(5.2) 



a,b 



where T is the kinetic-energy operator. Vat is the inter- 
action between the constituent particle a and 6, and the 
^Pauii is the Pauli projection operator given by Ea (|3.6p . 
The total wavefunction is described as a sum of ampli- 
tudes of the rearrangement channels (c = 1 9) of Fig. [3] 



c=ln,N,iyl,L,\S,I,K 

X 5„[$(ai)$(a2)[xi(")Xi(S')]5 
X [kf(rc)V^^^(Rc)],e^Pc)] 



(c) 

nlNLvXSIK 



JM 



(5.4) 



Here the operator Aj\f stands for antisymmetrization be- 
tween the two neutrons. 5^, xi('^~) ^'^d 7yi(S~) are 
defined already in Sec. III. 

The Pauli principle involving nucleons belonging to a 
and x{— 71, a) is taken into account by the orthogonality 
condition model (OCM) The forbidden states in 

Eq.dSH) are / = OS' for n and / = OS, IS, for x = a. 

We employ the VaN potential given in Ref.fl^l and 
the AV8 potential [2^ for the two-neutron parts. The 
ann [aan] binding energy derived from these potentials 
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(a)!-H (ESC) 



(b)2-H (ND) 



MeV 
-0.96 



H(aS-)cai+n+n,'-' 



-1.71 



1/2*- 



-3.06 



-0.57 .-■ 



-0.19 



-a+n+n+£. 



6He+S- 



-1.29 



-1 .83 



MeV 



sH(aH-)cai+n+n 



-0.57 



-0.96 



1/2*- 



-2.52 



-0.32 



-0.15 



-1.50 



-2.02 



-Ot+n+n+i 



6He+S- 



(i) kp =0.9 



I kp =1.055 (iii)kF=1.3 



(i)kF=0.9 (ii)kF =1.025 (iii)kF=1.3 



FIG. 4: (a) Calculated energy levels of -LH for three fcp values using ESC. (b) Calculated energy levels of ^-H for for three 
fcp values using ND. The energies are shown when the imaginary part of the dST interaction is switched off. The energies are 
measured from the a -I- n -I- n -I- H~ breakup threshold. The dashed lines are threshold. 



is less(over) bound by about 0.3 MeV (1 MeV) in com- 
parison with the observed value. Then, in calculations 
of the anriEr and aanEr four-body model, the central 
part of Van is adjusted so as to reproduce the observed 
ground state of ^He and ^Be. The Vaa and Vqb.- are the 
same as those in aafEr four-body calculations. As for 
the S^n parts, we employ the simple three-range Gaus- 
sian potentials derived from ESC and ND. The details of 
these potentials were already mentioned in Sec. II. Thus, 
in our treatments of anriEr and aavEr four-body sys- 
tems, ground-state energies of all subsystems of ann and 
OLOin are reproduced well. 



B. Results for g^H (annH ) 

Here we describe the results of the four-body calcula- 
tions for 2^H(annS-) with (T, J'^) = (3/2,1/2+). The 
basic question is whether this state is bound or not: The 
^He core is composed of an a and two weakly-bound 
('halo') neutrons. Due to the weakness of the S~n in- 
teraction, the binding between ^He and S~ is to a large 
extent determined by the dBr interaction. 

The calculated energies in the 1/2+ ground state are 
demonstrated in Fig. 0] as a function of fcp, for the two 
EN potential models without the imaginary part of the 
dBr interaction. These 1/2+ states are composed of 
the ground-state 0+ configuration of ^He coupled with 
the Os-state E" particle. The Coulomb interactions be- 
tween a and are taken into account. In the figure, the 
dashed lines show the positions of threshold energies of 
a-|-n-|-n-|-S~, ^He-|-5~ and _^H(aS~)cai-l-'^-l-n, respec- 
tively. One should be aware that the _^H(aS^)cai-l-"--l-fT- 
threshold energy depends on the fcp value of the adopted 
EN interactions. This situation is unavoidable, because 
the calculated energies for have to be used instead 
of the unknown experimental value. We see that in the 
case (i) fcp = 0.9 fm~^ with ESC the lowest thresh- 



old is _'^H(aS^)cai + n + n, and in the other cases the 
^He + S^ threshold is lower than the -5. H(aS~)cai + + 
threshold. On the other hand, in all fcp cases with ND 
the lowest threshold is ^Hc + E^ . The order of the 
_5.H(aS~)cai + n + n and ®Hc + threshold is deter- 
mined by the competition between a-E~ correlation and 
the a -{nn) correlation. 

More detailed results are given in Table IVIIIl where 
the calculated values of the conversion widths F and the 
aE~ and an r.m.s. radii are also listed. 

Now, let us compare the results for ESC and ND in 
the cases in which the imaginary part of the aE^ inter- 
action is switched off. As found in Table IVII (values in 
parentheses), the obtained aE~ states for ESC are more 
bound than those for ND (-1.71 MeV vs. -0.57 MeV 
for kp = 0.9 fm~^). In the annE~ system, however, the 
energy difference between ESC and ND becomes small in 
comparison with that in the aE~ system (—3.06 MeV vs. 
-2.52 MeV for fcp = 0.9 fm^^), as shown in Fig. Hand 
Table IVIIIl (values in parentheses) . This is because the 
'^^S'o and "^^Si nE" interactions of ND are more attrac- 
tive than those of ESC, as shown in Fig. [1] The stronger 
nE~ attraction in ND has the effect of a larger reduc- 
tion of the value of f^-'s- when one goes from the aE~ 
system to the annE~ system. 

Let us discuss the structure of _I H. In the case of fcp = 
0.9 fm-i for ESC, the lowest threshold is jLYl{aE-)^^y + 
n + n. Then, the particle is bound to the a particle 
mostly in the Os orbit, and the two valence neutrons are 
coupled to the aE~ subsystem. In fact, as shown in 
Table Ivnll f^s- is shorter than fa-n in this case. In 
other cases, the two valence neutrons are bound to the 
a core, and the E!~ particle is coupled to the ann(^IIe) 
subsystem, corresponding to where the r„_H- values are 
larger than the fa-n values. 

In order to see the structure of the _I H system visually, 
we draw the density distributions of (solid curves) 
and valence neutrons (dashed curves) in Fig. EJa) and 



10 



(a)^_H(ESC) 



10^ 

10-' 

10^ 

10-= 



a 

V 


(i) kp =0.9 


' ''\/ 






Xx 

1 1 1 X^i 1 1 



5 10 
r(fm) 



10^ 
10-' 

i. 1°"' 

^10-= 
10^ 

10-= 





(ii) kp =1.055 


/ ' Xn 
' ' Xv 






V 


n 

\\ / 






■^NL / 
/ 

^ ^x / 
^x / 

N \> 

X \, 

X \. 

\ \^ 






X \^ 
X \^ 

^. , , 





(Okp =0.9 
























\X^^ 




vx ^ 




\\ / 
\ \ / 




X / 

/ 








\ \ 

X \ 




X \ 
X 



15 5 10 

r(fm) 

(b)^_H(ND) 

(ii) kp =1.025 



15 



5 10 
r(fm) 



15 




(iii) kF=1.30 






(iii) kF=1.3 


/ / '^x\ 


n 


x\ 

; X 


X 

X 

X 

X 

X 

X 

, , 



15 



5 10 
r(fm) 



15 



FIG. 5: (a) Calculated density distribution of q, H~ and a valence neutron for three ftp values using ESC. (b)Calculated density 
distribution of q, H~ and a valence neutron for three fcp values using ND. The wavefunctions of aST without the imaginary 
part of the dST interaction are used. 



(b) for ESC and ND, respectively. For comparison, also 
a single-nucleon density in the a core is shown by a dot- 
ted curve in each case. It turns out, here, that as the 
binding energies of -IH become smaller, the density 
distribution has a longer tail. As is well known, ^He is 
a neutron-halo nucleus. It is interesting here to see the 
overlapping of the S~ distribution with the halo-neutron 
distribution. In the case of fcp = 0.9 fm^^ for ESC, since 
the lowest threshold is _5.H(aS~)cai -I- n -I- n, the density 
of the S~ particle has a shorter-ranged tail than that of 
the two valence neutrons, but is extended significantly 
away form the a core. This situation can be visualized 
as three layers of matter distribution, the a core, a "Er 
skin, and neutron halo. When the lowest breakup thresh- 
old is ^He+S^ , the S~ density is longer-ranged than that 
of the valence neutrons due to the weaker binding of the 
S~ particle. Then, the density distribution of _IH shows 
the three layers of the a core, neutron halo, and halo. 
Namely, a double-halo structure of neutrons and ex- 



ists, in which the attractive Coulomb interaction plays an 
essential role. These features can be considered as new 
forms in baryon many-body systems. 



Table I Villi lists the binding energies of the _I H sys- 
tem calculated with and without the Coulomb interac- 
tion for each fcp value. For ESC (fcp — 0.9 fm^^) and 
ND (/CF = 0.9 and 1.025 fm^i) the ground states of H 
are found to be weakly bound states, when the Coulomb 
interactions are switched off. Therefore, the -IH sys- 
tems are seen to have nuclear-bound states, if we take 
reasonable values /cp < 1 fm^^. This means that an ex- 
perimental finding of a -IH bound state indicates the 
existence of an a'BT bound state in which the even-state 
spin-independent part of the SA^ interaction is substan- 
tially attractive. This statement is almost independent 
of the interaction model. 
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(a) ;°Li (ESC) (b) H-Li (ND) 

OMeV _ OMeV 

a+a+n+ E a+a+n+ ^ 



rLQ4 

-2.89 -2.92 



9Be.- ^^^-9Be. 

-1.57 -1.57 _■! 83 



"^■^^ =aT3 |Li(aaSi,„+n-2.89 -2.97 



? Li (oa S -),3|+n ,/' , zlZL-;. 



^.43 ,-' 



5.17 // -5.16 
7.83 //' 



2- -5.85 . 
■I' -5.99 



-4.53 



r 

2 -8.35 

(i)kp=0.9 (ii)kp =1.055 (iii)kp=1.3 {i)kF=0.9 (ii)kp =1.025 {iii)kF=1.3 

FIG. 6: (a)Calculated energy levels of i-Li for three fcp values using ESC. (b) Calculated energy levels of g^Li for three fcp 
values using ND. The energies are shown when the imaginary part of the qH~ interaction is switched off. The energies are 
measured from the a + a + n + S~ threshold. The dashed lines indicate thresholds. 



C. Results for Li (aan^ ) 

The calculated results for i!?Li(aQ;nS~) within the 
four-body model are displayed in Fig. [6] for the lowest 
T = 1 doublet state energies (J'" = 2", 1"). The 3/2" 
ground state of the core nucleus ^Be is bound by about 
1.57 MeV with respect to the a + a + 7i threshold. We 
emphasize that, if the aaS" system is bound as shown in 
SecJV, then the i°Li(aanE!~) system is surely expected 
to be bound because the interaction between the and 
a p-orbit neutron is weakly attractive. 

Although the binding energies of i°Li (J^ = 2^, 1^) 
are found to be fairly sensitive to the choice of the fcp val- 
ues, especially, in the case of ESC, we think the results 
with kp ~ 1.0 fm~^ are most acceptable. It is interesting 
to note that the _^Li(aQ;S^)cai + n threshold comes be- 
low the ^Be-|-S~ threshold in most cases. It is reasonable 
that the lowest breakup threshold is _^Li(aaS^)cai + n, 
because the value oi kp in the A = 10 system has to 
be ~ 1.0 fm~^, similar to the i?Be system. It is no- 
table that the binding energies of Li measured from the 
_?_Li(Q:Q;S~)cai + n thresholds are similar to each other 
hi both cases of ESC(2.6 MeV) and ND (2.7 MeV). 

We expect such structure, a valence neutron cou- 
pled to -^Li hypernucleus, since the lowest threshold is 
-?.Li(aQ;S~)cai + n. While, if the fcp value becomes by 
chance much larger, then the S" particle is coupled to 
the ground state of ^Be. Because the lowest threshold is 
^Be -I- S~ [See case (iii) in Fig. [5]. 

The 2" (1") state is dominated by the ^^Si {^'^Sq) 
component of the two-body nS~ interaction. As demon- 
strated in Fig. [H the ^^Si interaction for ESC (ND) is 
more (less) attractive than the ^^5*0 interaction. There- 
fore, the 2~ (1~) state of i° Li becomes the ground state 
in the case of ESC (ND), as shown in Fig. [6] 

More detailed results are given in Table HXl which lists 
also the calculated values of the conversion widths F and 



the r.m.s. radii, r^..^- and rQ,_„. We show here the 
results with and without the aS~ Coulomb interaction. 

As seen in Table IIXI the decay widths T calculated 
with ESC are much larger than those for ND, mainly 
because the ^^5*0 SA^-AA coupling interaction in ESC is 
far stronger than that in ND. 

The r.m.s. distance, fa-s- , both for ESC and ND, are 
comparable to the fa^„ values in the cases of choosing 
plausible kp values. In order to illustrate this situation 
visually, we show the density distributions of (solid 
curves) and a valence neutron (dashed curves) in Fig. [7] 
(a) and (b), where a single- nucleon density in the a core 
is indicated by dotted curves. The upper part (a) and 
the lower part (b) in the figure are for ESC and ND, re- 
spectively. In the {K~,K^) reaction, if the spin-nonflip 
transition dominates, the 2~ state of i°Li is selectively 
excited. Then, we show the density distributions of 2~ 
states. The densities of the S~ and a valence neutron are 
extended significantly away from the a core. 

Let us see the effect of the Coulomb interaction be- 
tween the a and . In Table 11X1 we list the binding 
energy of the 2~ state in i°Li with and without the 
Coulomb interaction. The most important result is that 
the states obtained for i° Li are bound without Coulomb 
interactions, namely as nuclear-bound states, for all k-p 
values for both ESC and ND. 

According to the above calculations, we can surely ex- 
pect the existence of the nuclear bound state with the 
predicted S~ binding energies of -B^- = 3.26 MeV (ESC, 
kp = 1.055 fm-i) and 2.85 MeV (ND, kp = 1.025fm-i), 
when the imarginary part of aS~ interaction is taken ac- 
count of. These i?^- values seem to be a little smaller 
in comparison with the empirical value ~ 4.5 MeV sug- 
gested in the ^"^ClK' , K~^)^'^Be reaction. However, to 
produce i!?Li, we propose to perform the ^'^B{K~ , K^) 
reaction experiment at J- PARC in addition to that with 
a ^^C target. 



12 




FIG. 7: (a)Calculatcd density distribution of a, H~, and a valence neutron for three fep values using ESC. (b)Calculated density 
distribution of q, "BT , and a valence neutron for three fcp values using ND. The wavefunctions of aH~ without the imaginary 
part of the dST interaction are used. 



We say that the aanEr (i" Li) system produced by the 
(K~ ,K^) reaction on is suitable to investigate dBr 
interactions, namely the spin-independent terms of even 
and odd-state SAT interactions. 



VI. SUMMARY AND OUTLOOK 

In anticipation of priority experiments to be done at 
the J-PARC facility, we have carried out detailed struc- 
ture calculations for several light p-shell S-hypernuclei, 
i?Be, -'iH, -'^Li, -IH and i°Li, in order to investi- 
gate whether we can expect the existence of bound states 
of the S~ hyperon. The calculational framework is mi- 
croscopic three- and four-body cluster models using the 
Gaussian Expansion Method which has been proved to 
work quite successfully in obtaining reliable numerical 
solutions. 

One of the essential issues in preparing such detailed 



calculations is what kind of SA^ interactions one should 
use, because there are no definitive experimental data 
for any S-hypernucleus, and also because there are large 
uncertainties in the spin and isospin dependence in the 
existing SA'' interaction models. The only existing ex- 
perimental indication, from the ^^C(i4r~, i<r"'")J?Be re- 
action spectrum, is that the ^^B-S" interaction is sub- 
stantially attractive. However this constraint is helpful 
in excluding most of the iSJ/a-invariant BB interaction 
models which lead to repulsive S- nucleus potentials. In 
this work, we used two 'E.N potential models, ND and 
ESC, which give rise to substantially attractive S-nucleus 
potentials in accordance with the experimental informa- 
tion. Although the spin- and isospin-components of these 
two models are very different from each other due to 
the different meson contributions, we can reliably speak 
about the spin- and isospin-averaged properties such as 
(5(±) = (g[,J' + 3G^^) + 3G1^) + 9G1^))/16. This is 
why we have focused our attention on the a-cluster based 
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TABLE Vlll: The calculated binding energies of 1/2+, E and 
r.m.s radii, f^_-^- and fa-n , in the 2-H(annH~) system for 
several values of kp. The values in parentheses are energies 
when the imaginary part of the aH^ interactions are switched 
off. The energies are measured from the a + a + n + E~ 
threshold. 

(a) glHCESC) 



fcF(fm-i) 0.9 



with 

Coulomb E (MeV) 
r(MeV) 
fa_„(fm) 



1.3 



1.055 

-2.76 -1.63 -1.22 

(-3.06) (-1.83) (-1.29) 

2.64 1.15 0.31 

3.68 5.58 9.92 

4.04 4.11 4.19 



without E (MeV) 
Coulomb 

r(MeV) 



— 1.68 unbound unbound 
(-1.96) (-1.09) (unbound) 
2.09 



(b) slH(ND) 



with fcF(fm ) 

Coulomb E (MeV) 

r(MeV) 

fg-n (fm) 



0.9 1.025 
-2.51 -2.01 
(-2.52) (-2.02) 
0.27 0.15 
4.48 5.35 
3.92 3.99 



1.3 
-1.50 
(-1.50) 
0.032 
7.55 
4.11 



without E (MeV) -1.62 -1.26 unbound 
Coulomb (-1.63) (-1.26) (unbound) 

r(MeV) 0.22 0.10 



systems and started with an investigation of the nuclear 

spin- and isospin-saturated systems such as aS~(_^H) 
and aQ;S~(„?.Li), so as to get a firm basis of our analy- 
ses. 

However, the pure a-clustcr systems such as aS" 
and aa'E^ cannot be produced directly, because there 
are no available nuclear targets for the {K~,K^) re- 
action. Thus, in order to explore realistic experimen- 
tal possibilities, we have extended the calculation to the 
four-body !B~-systems having one or two additional neu- 
trons. This explains why we took the _lH(annS~) and 
~-Li{aanE~) hypernuclei as the typical S~-systems in 
this paper. 

The major conclusions are summarized as follows: 
(1) In order to be consistent with the existing exper- 
imental indication that the S-nucleus interaction is at- 
tractive, the fine tuning of the ND and ESC potential 
models has been made for applications to S-hypernuclei 
by adjusting the hard-core radius Vc in ND and the ay 
parameter for the medium-induced effect in ESC, respec- 
tively. Then the EN G-matrices were derived and rep- 
resented in terms of three-range Gaussians with the kp 
parameter expressing its density-dependence within the 



TABLE IX: The calculated binding energies, E of the 1]~ and 

2]~ states in the i° Li(aQnE~) system for several values of fcp. 
The values in parentheses are energies when the imaginary 
part of the aH^ interactions arc switched off. The energies are 
measured from the a + a + n + E^ threshold. The calculated 
r.m.s. radii, Va-s- 1 ^a-n and fa-a of 2~ state using ESC 
and ND. 





(a) i°- 


.Li(ESC) 






with Coulomb 


A;F(fm-i) 


0.9 


1.055 


1.30 


2" 


E (MeV) 


-7.99 


-4.83 


-2.87 






(-8.35) ( 


-5.16) ( 


-3.13) 




r(MeV) 


5.87 


3.63 


1.71 




^a-B- (fm) 


3.05 


3.72 


5.03 




ra_„(fm) 


3.55 


3.70 


3.83 




fa-a(fm) 


3.25 


3.41 


3.54 


1" 


E (MeV) 


-7.48 


-4.42 


-2.64 






(-7.84) ( 


-4.77) (■ 


-2.89) 




r(MeV) 


5.72 


3.44 


1.50 


without Coulomb 


E (MeV) 


-5.54 


-2.76 


-1.41 


2" 




(-5.93) ( 


-3.14) ( 


-1.63) 




r(MeV) 


5.39 


3.00 


1.10 




(b)i^ 


LLi(ND) 






with Coulomb 


/cF(fm-') 


0.9 


1.025 


1.3 


2" 


E (MeV) 


-5.83 


-4.42 


-2.92 






(-5.85) ( 


-4.43) ( 


-2.92) 




r(MeV) 


0.75 


0.42 


0.10 




f„_H- (fm) 


3.55 


4.10 


5.40 






3.64 


3.72 


3.83 




fa-Q(fm) 


3.35 


3.44 


3.54 


1" 


E (MeV) 


-5.98 


-4.53 


-2.97 






(-5.99) ( 


-4.53) ( 


-2.97) 




r(MeV) 


0.77 


0.43 


0.10 


without Coulomb 


E (MeV) 


-3.75 


-2.60 


-1.54 


2" 




(-3.76) ( 


-2.61) ( 


-1.54) 




r(MeV) 


0.62 


0.32 


0.005 



nucleus. The AA — EN — SS coupling term in the bare 
interaction is renormalized into the imaginary part in the 

G-matrix interactions 

(2) First we performed the aatE~ four-body calcula- 
tion of J?Be and found that kp = 1.055fm~^ for Gpsc 
and 1.025 fm"'^ for Gnd arc most appropriate to pro- 
duce the B2-^{J = 1") = 2.2 MeV (without Coulomb 
interaction) which is that suggested empirically. These 
values around kp ^ 1-0 fm"-'^ are reasonable, because 
they agree roughly with the values estimated from the 
average density in a A = 12 nucleus. Then, we naturally 
allow smaller kp values for smaller mass numbers. 

(3) In the basic structure calculations for aS~(g.5.H) 
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and aaE~{JLLij systems, we have tested three vahies 
of kp parameters, kp — 0.9,1.055 and 1.3 for ESC and 
kp = 0.9, 1.025 and 1.3 for ND, respectively. In the 
aS~ system, for which kp ~ 0.9 fm~^ is considered to 
be reasonable, we obtained only Coulomb-assisted bound 
states with small binding energies, since they disappear 
without the Coulomb interaction. In the aaS" system, 
on the other hand, nuclear bound states are obtained for 
the acceptable range of kp between 0.9 fm~^ and 1.05 
fm~^. The calculated binding energies of ESC are larger 
than those of ND, and also the fcp-dependence is more 
sensitive in ESC. If these predictions are confirmed, di- 
rectly or indirectly, in future experiments, then it will 
provide a good check for the spin- and isospin-averaged 
SiV interaction strengths. 

(4) For the lightest realistic example, -lll(annS^), 
the four-body calculation predicts the existence of nu- 
clear bound states in both cases of ESC and ND at rea- 
sonable kp values of around 0.9 fm^^. It is interesting 
to note that the addition of two neutrons to the a5^ 
system gives rise to about 1.3 (2.0) MeV more binding 
for the ESC (ND) cases, respectively. If the experiment 
is carried out to observe the -IH bound states, it is use- 
ful to extract information about the even-state spin- and 
isospin-averaged part of the "BN interaction acting be- 
tween the a and S~. 

(5) For the second realistic example, the i°Li(aanS~) 
hypernuclus, we have obtained the nuclear bound 
states as a result of careful four-body calculations with 
kp ~ 1.0 fm~^. This result is essentially based on the 
averaged attractive nature of the SiV interactions acting 



in the three-body subsystem of aa:i . It is remarkable 
to have similar binding energies of the J — 1~ state for 
both the ESC and ND interactions (-4.8 MeV vs. -4.4 
MeV with respect to the a -I- a + n + threshold). The 
order of the doublet states (J = 1^, 2^) is calculated to 
be opposite for ESC and ND as a result of the difference 
in the '^^S\ and ^^Sq components of ESC and ND acting 
between the S~ and a neutron. 

In conclusion, it will be quite interesting to observe 
the newly predicted bound states in future {K^ ,K^) 
experiments using the ''Li and -'^'^B targets in addition to 
the standard ^^C target. Experimental confirmation of 
these states will surely provide us with definite informa- 
tion on the spin- and isospin-averaged SiV interactions; 
note the information on its even-state part from oiEr and 
adEr and its odd-state part from aaS~. Such a plan 
is a challenging project in the study of S-hypernuclei 
that have yet to be explored. In order to convert the 
present predictions into concrete experimental proposals 
at J-PARC, the reaction cross sections should be esti- 
mated for the ^Li(is:-,X+)^^H, iOLi(if-, if+)io Li and 
i2C(if - , K"-)^- Be reactions. 
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